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The unsteady, sonic, small disturbance potential equation 
was approximated by reduction to parabolic form and solved by 
Laplace and Fourier Transform and Collocation techniques, 
including thickness effects through the Oswatitsch parameter. 
The Transform techniques produced reliable solutions and the 
Collocation solution was usable with qualifications, for some 
parameter values, stemming from the choice of approximating 
series. Stability boundaries for an unstaggered cascade 
oscillating in sonic flow were calculated with the Fourier 
Transform technique revealing that blade interference is 
strongly destabilizing. The Collocation technique was suc- 
cessfully applied to the two-dimensional, unsteady, porous 
walled, sonic wind tunnel interference problem, heretofore 
unapproached with the inclusion of thickness effects. 
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LIST OF SYMBOLS 



K 2 /(X 2 + 4K 2 ) 

nondimensional position of pivot location, 
a divided by reference length 

position of pivot in pitch oscillation 

variable in Fourier Transform solution, functions 
of interblade phase angle or can be = XP^/4 

nondimensional pivot location, b divided by 
reference length 

airfoil or blade semichord 

dimensionless speed of sound; speed of sound 
divided by reference velocity 

lift coefficient, positive upward 

moment coefficient, positive leading edge up 

pressure coefficient 

airfoil or blade chord 

variable in Fourier Transform solution, function 
of interblade phase angle 

a 

e 

KP 2 /4 

constant in solution of Fourier transformed 
potential equation 

constant in solution of Fourier transformed 
potential equation 

first integral of Y, used in Fourier Transform 
solution 

nth integrals of V, n > 1 

variable necessary in Fourier Transform 
solution = 1,2 
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H(X,T) 

H(Z) 

h 

h 

h(X) 

h 0 

h i 

*n 

I(X) 

i 

J(X) 

K 

K 1 

L 

M 

n 

P 

P 

Q 

R.P. 



S 



N 



nondimensional Y position of airfoil or blade 
surface 

Heaviside step function 

nondimensional interblade distance, P/C 

nondimensional mean Y position of airfoil or 
blade surface 

nondimensional amplitude of oscillatory airfoil 
or blade surface 

amplitude of plunge oscillation 

nondimensional amplitude of plunge oscillations, 
divided by reference length, i = 1 or 2 

integrals appearing in Laplace Transform and 
Collocation solution 

integral values for Fourier Transform solution 
imaginary constant, or variable index 

integral values for Fourier Transform solution 
reduced frequency; wb/U^ • 
u)c/U^ 

Oswatitsch parameter, X, or number of collocation 
points to be used 

Mach number of reference velocity 
variable index 

interblade distance, or variable in I n 
Laplace Transform transformation variable 
X + 2iK 
real part 

variable of integration for Laplace Transform and 
Collocation solutions 

see T, t 
N 
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T = dimensionless time, t multiplied by reference 

velocity and divided by reference length 

T m = integral values necessary in Fourier Transform 
solution 

t = time 



U = X-S, variable used in evaluation of I 

n 

= reference velocity, freestream or critical 



u m 



interference source distribution, m = 0 or 1 for 
reference or adjacent blade, respectively 



v = nondimensional downwash velocity at blade chordline, 
used in Laplace Transform and Collocation solutions 

v m = source distribution on airfoil or blade, m = 0,1, 
for reference or adjacent blade, respectively 

w(X) = nondimensional downwash velocity at blade chordline, 
used in Fourier Transform solution 



X = dimensionless coordinate in freestream direction; 

x divided by reference length 

Xq = nondimensional pivot location, a divided by b 

x 1 = X + 1 

XA = upper limit of 'correction integral' 
x = coordinate in freestream direction 



Y = dimensionless coordinate perpendicular to stream 

direction 



y = coordinate perpendicular to stream direction 

cxq = amplitude of pitch oscillation 
Y = ratio of specific heats 

= perpendicular coordinate from blades, i=0,l 
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variable in summation for <J(X) 

coefficients of isolated blade source distribution 
series, m defines blade (reference or adjacent), 
n defines sequence number in series 

Oswatitsch parameter, (y+l)4i^. 

variable in transformed potential equations 

Fourier Transform transformation variable 

variable in summation for I(X) 

coefficients of interference source distribution 

series, see 0 

5 mn 

coordinate along airfoil or blade chord from 
midchord or leading edge, i = 1 or 2, respectively 

nondimensional interblade distance, P divided by b 

interblade phase angle 

porosity parameter 

airfoil or blade thickness/chord ratio 

dimensionless unsteady two-dimensional velocity 
potential, potential divided by reference length 
and reference velocity 

nondimensional perturbation potential, or dummy 
variable 

isolated blade potential for the reference blade 
isolated blade potential for the adjacent blade 
interference potential for the reference blade 
interference potential for the adjacent blade 
steady flow perturbation potential, nondimensional 
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i| ; = nonsteady flow perturbation potential, nondimensional 

ijj = transformed nonsteady perturbation potential 

to = circular frequency of pitch or plunge oscillation 
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I. INTRODUCTION 



Nonsteady, interfering, transonic flows are of increasing 
importance to aeroelasticians . Recent design concepts for 
transonic transports, high maneuverability fighters, biplanar 
space shuttle configurations, and so called transonic axial 
flow compressors, to name a few, involve such flows when aero- 
elastic problems are considered. Ironically, while sub/super- 
sonic flutter problems have been extensively studied, the 
transonic problem has been circumvented by either rapid 
acceleration through the regime or operation avoiding it. It 
is now apparent that the question of the effect of transonic 
flow, including possible adjacent flov; field interference, 
on flutter must be confronted in a practical way. 

Of course, even without interference, nonsteady transonic 
flows challenge the aeroelastician concerned with aerodynamic 
damping of structural vibrations. In that regime, pressure 
perturbations generated by a body and propagated at the local 
speed of sound tend to accumulate near the body itself. Thus, 
the magnitude of aerodynamic forces is increased and, perhaps 
more importantly, a large phase lag can exist between the 
forces and the motion producing them. It is this phase lag 
that can lead to negative aerodynamic damping and, potentially, 
dangerous flutter. 

It is not enough to analytically study the flutter charac- 
ter of particular configurations. Wind tunnel experiments 
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must be conducted to verify analytical results before a 
configuration can be committed to flight test. However, 
flow perturbations from a wind tunnel test model may reflect 
from the tunnel walls back onto the model requiring a correc- 
tion to test results due to interference. In steady transonic 
flow, such wall interference has been thoroughly examined and 
porous or slotted wall tunnels have been utilized for some 
time to reduce its effect. In nonsteady flow, the porous 
wall boundary condition is a difficult one for classical 
analytical approaches to handle. Only Drake [Ref. 6] has in- 
cluded wall porosity in a theoretical analysis of a thin air- 
foil oscillating in a transonic wind tunnel and his results 
are by no means complete. Thus, there is a great need for 
solutions that can yield the magnitude of unsteady wall 
interference and the effect of wall porosity upon it. 

Blade flutter in transonic axial compressors is a complex 
problem involving, among other things, interference of highly 
perturbed transonic flow fields. Historically, sub and super- 
sonic axial compressor flows have been simulated by a series 
of airfoils (or blades) placed in close proximity (a cascade) 
so that the two dimensional flow models that in a rotating 
compressor stage. For the transonic compressor stage the 
mixed flow character of the real flow makes the two-dimen- 
sional approximation questionable, e.g., McCune [Ref. 1 6 ] . 

It must be recognized that accurate modeling of such a complex 
flow field is well beyond the state of the art. However, the 
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cascade approximation may yield important information relative 
to the trends of blade flutter. 

Nonsteady transonic flow analysis begins with the deriva- 
tion of the nonlinear small perturbation potential equation 
which Lin, Reissner, and Tsien [Ref. 15] and Landahl [Ref. 

13] have shown can be linearized for sufficiently high fre- 
quency. Since aeroelastic analyses are, generally, concerned 
with the stability of small amplitude oscillations, the solu- 
tions to the nonlinear potential equations shall be written 
as the sum of a steady component and a much smaller nonsteady 
component. The result is two equations, the usual nonlinear 
steady flow equation and a linear nonsteady flow equation 
containing variable coefficients depending upon the steady 
flow . 

Transonic cascade flow for thin blades and in-phase oscil- 
lations was approached by Hamamoto [Ref. 10] using Fourier 
Transform techniques on the linearized potential equation. 
Gorelov [Ref. 8] solved the linear nonsteady equation for 
mildly supersonic flow past thin blades with a superposition 
of isolated blade and interference potentials represented as 
trigonometric series. See also Samoylovich [Ref. 22]. 

Elder [Ref. 7] revised Gorelov’s Collocation technique by 
using a power series and solved for the steady sonic flow 
cascade potential and the coefficients of the series for the 
nonsteady flow, out-of-phase oscillation case. Elder made use 
of the Oswatitsch [Ref. 20] parameter (used by Oswatitsch for 
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steady flows), as extended by Teipel [Ref. 30] and Hosokawa 
[Ref. 11] for nonsteady flows, to allow for blades of nonzero 
thickness (nonuniform steady potential flow). 

Sonic wind tunnel interference in steady flow was solved 
by Sandeman [Ref. 23] using Spreiter's local linearization 
technique, this after Guderley [Ref. 9] and Spreiter et al 
[Ref. 27] showed that the choked tunnel results can be cor- 
rected so as to simulate unbounded flows. Murasaki [Ref. 1 8 ] 
solved steady transonic tunnel flows for porous or slotted 
walls. Miles [Ref. 17] and Drake [Ref. 5] solved the super- 
sonic nonsteady solid and porous wall cases, respectively. 
Drake [Ref. 6], then, utilized Laplace Transform techniques 
and the linear potential equation for a thin airfoil between 
porous walls for transonic flows. Recently, Savkar [Ref. 2^4] 
has Laplaced transformed the linearized equation and studied 
airfoil stability in a solid walled wind tunnel. 

Wind tunnel interference can be approached, according to 
Sears [Ref. 25], with suitably located sensors within the 
test section (not on the model) and variable wall shape or 
porosity and/or plenum pressure behind the wall. The sensors 
would modify the tunnel parameters until unbounded flow con- 
ditions were present at them. Sears argued that unbounded 
flow conditions are thus obtainable upon the model. 

This work began with the object of applying a Collocation 
technique to the porous wall wind tunnel problem with the con- 
viction that that method could suitably handle the classically 
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difficult porous wall boundary condition. While proceeding 
toward that goal, an efficient approximate solution was de- 
veloped for studying the aerodynamic stability of the sonic 
oscillating cascade and stability boundaries were calculated 
for the cascade as well as for a biplanar configuration. 

To accomplish the initial objective of this work, the 
nonsteady potential equation was reduced with the aforemen- 
tioned Oswatitsch parameter to approximate the variable 
coefficient therein. The resulting reduced equation was 
solved by a Collocation technique for the cascade problem 
following Elder [Ref. 7]. For numerical checks upon the 
Collocation technique, the equation was Laplace transformed 
extending the steady flow work of Sandeman [Ref. 23 ] and 
Fourier transformed in a manner analogous to Hamamato [Ref. 
10], These extensions were made flexible enough to handle 
blade thickness variations (through the Oswatitsch parameter) 
and variable interblade phase angle. 

The Elder power series initially utilized in the Collo- 
cation solutions to the cascade problem was found to have 
numerical difficulties that increased computation time re- 
quirements and reduced accuracy. Thus, a 'step-wise-linear' 
series was introduced and the resulting Collocation method 
solution to the cascade problem agreed well enough with the 
other two solutions to allow some confidence in its subsequent 
adaptation to the porous wall wind tunnel problem. 
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The Collocation solution was found to be quite flexible. 
It can be altered to allow streamwise variation of the Os- 
watitsch parameter, in a manner reminiscent of Spreiter's 
local linearization, to improve accuracy. Further, it allows 
a variation in wall porosity for possible simulation of the 
technique mentioned by Sears. 
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II. UNSTEADY TRANSONIC FLOW THEORY 



A. EQUATIONS FOR GENERAL OSCILLATORY MOTION 

The governing equation considered below arises from the 
assumption that the flows of interest are of frictionless 
perfect gas. Thence, the equation of state and Eulerian 
equations of motion are invoked and simplified, assuming the 
absence of heat conduction and strong shock processes. Fur- 
ther, a velocity potential is introduced which reduces the 
number of dependent variables at the imposition of the re- 
quirement of flow irrotationality . Finally, a small distur- 
bance analysis, including the unique properties of transonic 
flows (see Refs. 12 and 13) , yields, in two dimensions 



This equation is applicable, approximately, to flows containing 
weak shocks, incidentally (see Ref. 13). 

In the study of unsteady aerodynamics, it is customary 
and convenient to expand the solution into steady and 
oscillatory nonsteady components. Thus, 



4> tt + 2< J ) xT + ^ ^2 




4>(X,Y,T) = 4>(X,Y) + R.P. !>(X,Y)e lKT ], 



( II-2 ] 



where K is the reduced frequency. 
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Substitution in equation (II-l) yields a coupled partial 
differential equation, 

(1-M 2 H xx + Jp YY - M 2 (y+1)[^ x ^ xx + ^ x ^ xx 

+ M 2 K 2 i|j - 2iM 2 K^ x = 0. (II-3) 

This equation is nonlinear, of mixed elliptic-hyperbolic 
type, and contains variable coefficients. However, because 
flutter analysis is, generally, concerned with the stability 
of small perturbations about a steady flow, it is appropriate 
to assume that the oscillatory nonsteady component is small 
compared to the steady flow potential. Thus, the governing 
equation becomes 

(i-m 2 H xx + ip yy - m 2 (y+i)U x i|, xx + ^ xx ^ x l 

- 2iM 2 K^ x + M 2 K 2 ^ = 0, (II- 1 *) 

with the usual nonlinear steady flow equation also resulting. 
Equation (11-4), though linear, retains variable coefficients 
dependent upon the steady flow solution, the nature of which 
(subsonic or supersonic) determines its type (elliptic or 
hyperbolic ) . 

B. REDUCTION OF THE GOVERNING EQUATION 

First, the reference Mach number shall be unity, reducing 
equation (II-4) to 
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^YY “ ^ Y + 1 ) C ‘J’x'^XX + ^XX^X^ ” + ^ 2 iJ> = 



( II— 5 ) 



Second, the approach taken by Teipel [Ref. 30] and 
Hosokawa [Ref. 11] is utilized. That is, the nonconstant 
coefficients on and \p^ are approximated as follows. 

First, the steady flow velocity perturbation is assumed to 
be small in the X direction, implying that the airfoil or 
blade is of small thickness, 

\ 1 0. (H-6) 

Second, the steady flow fluid acceleration over the airfoil 
or blade is assumed to be a constant (parabolic arc airfoils 
exhibit nearly this characteristic), which allows the intro- 
duction of the Oswatitsch parameter 



X 



(y+i)4>xx 



constant > 0. 



(H-7) 



The Oswatitsch parameter was originated by Oswatitsch and 
Keune [Ref. 20] for steady flottf, incidentally. Thus, the 
governing equation has, finally, become 

^ yY - (A+2iKH x + K 2 ^ = 0. (II-8) 

Note that this equation is of the parabolic type. 
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C. BOUNDARY CONDITIONS 



1 . Flow Tangency 



The usual surface boundary condition is 




(H-9) 



where F(X,Y,T) describes the surface of a body moving as a 
function of time. The boundary condition requires flow tan- 
gency , that is the total derivative of F is zero. Then, con- 
sidering a thin airfoil restricted to vanishingly small ampli- 
tude out of plane oscillations (recall the small nonsteady 
potential argument), its upper and lower surfaces can be 
described by 



F u (X,Y,T) = Y - H u (X,T) 



(11-10) 



F & (X,Y,T) = Y - H & (X,T), 



where 



H U (X,T) = Th u (X) + R.P.[h u (X)e lKT ] . 



(II-ll) 



With the definition of <J>(X,Y,T) from equation (II-2) 



and H(X,Y,T) from above, the linearized flow tangency 
condition becomes 



ipy(X,0±) = + iKh(X) . (II-12) 

2 . Zero Upstream Disturbance 

The governing equation has become parabolic as shown 



by equation 


(II-8), requiring, upstream of the leading edges. 




<KX,Y) = 0. (II-13) 
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III. FORMULATION OF SOLUTIONS 



A. SOLUTION THROUGH THE LAPLACE TRANSFORMATION 



1. Cascade Geometry and Boundary Conditions 




k 2b — — — *| 

p 

* n 0 

a —— " " - - 






1 



Fig. 1. Unstaggered Cascade. 

Consider a pair of the parabolic arc airfoils or 
blades in an infinite unstaggered cascade (row) immersed 
within a transonic flow field. The blade chord is 2b, inter- 
blade spacing is P, the lower blade is the reference blade 
(variables associated with it are subscripted, zero), and 
the upper is referred to as the adjacent blade (subscript, 
one). The coordinate reference is fixed to the reference 
blade midchord for consistency with Teipel [Ref. 30] with 
whose work comparisons were made. Pitch and plunge ocsilla- 
tions can be described, pitch axis at a, magnitude of pitch 
dp , of plunge h . The adjacent blade performs a harmonic 
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pitch/plunge oscillation arbitrarily phased with the oscilla- 
tion of the reference blade. 

Two n coordinates are specified, with rig and 
related such that 



"l = "o + P - 



(m-i) 



The following nondimensional variables are defined: 



x ■ TT ’ 



Y — 

. , m 

m b 



P = vT 



= 0 , 1 , 



( III-2 ) 



X rt = =■ 



a 

b * 



h n = 



_ 0 



Note that the linear nondimensionalizing parameter 
will be the blade semichord for all variables in this analysis. 
Thus, for example, the reduced frequency 



y _ Wb 
K " * 



(III-3) 
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From considerations in Section II, the boundary con- 



ditions can be expressed, for the reference and adjacent 
blades, as 



8h 

’J’y (X , Y Q = 0 ) = v 0 (X) = - 3 ^ + iKh u , (III-4) 



and 



^ Y (X,Y Q =p) = V;L (X) = e ia v 0 (X) , (HI-5) 

where a is the interblade phase angle. 

Then, for pitch and plunge of small magnitude 

v Q (X) = - cx 0 [l+iK(X-X 0 ] + iKh . (III- 6 ) 

For computations for the steady flow case, a 

p 

parabolic arc airfoil is assumed (Y = t( 1-X ), see Ref. 30) 
and 



Vq (X) = - 2tX. 



(III-7) 



2 . Solution of the Oscillatory Sonic Flow Equation 
The governing equation from Section II is 



iJ’yy “ + K V = °j 



'X 



(III- 8 ) 
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where 



Q = X + 21K. 


(III-9) 


For utilization of the Laplace Transformation 
a new streamwise coordinate 


define 


x 1 = X + 1. 


(III-10) 



Denoting the Laplace transformation of a function 
with respect to X-^ with an overbar, i.e.. 



- 'f *° " pX l 

f(p) = Z[f(X, )] - / e 1 f(X, ) dX n , 

x o - 1 x 

transform equation (II-8) to 


(III-ll) 


— 2— 

^ Y y - P ^ = 0 J 

where 


( III-12 ) 


2 2 
y = Qp - n . 

The boundary conditions become 


(III-13) 


= v Q (p) , 


(III-14) 
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(HI-15) 



= elO v 0 ( P) ' 



The solution to equation (III-12) subject to 
(III-14 and 15) is 



iKP,Y) = 



V 0 (P) e^ g cosh(yY) - cosh[y(Y-p)] 
y sinh( yp ) 



(III- 16 ) 



The hyperbolic functions are expanded, as follows, 
to allow a term by term inversion 



1 

sinh(yp) 



2 

(e pp - e" pp ) 



2e‘ pp Z e‘ 2ni,p 
n=0 



(III-17) 



then 



cosh( yY ) _ ”, e ~y[-Y+ (2n+l) p] + -y [Y+ ( 2n+l ) p ] , 

sinh(yp) „ =n l 



(HI- 18 ) 



and 



cosh[y (Y-p) ] 
sinh(yp) 



2 {e -y[2np+2p-Y] + e -y[2np+Y] } 
n=0 



(III-19) 



Thus , 



♦ (p,Y) - t: t {v Q (p)e io [e -p ^ -Y+2np+p) + e -vW+2np+p)-| 



p n=0 

-^ 0 (p)[e" ll( - Y+2np+2p) + e-W+Y+SnpJjj^ 



( III-20 ) 
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For inversion, recall the formula 



~ -p x i a/1(X i -jp 

J — dX 1 = -t/tT , R.P.[p] > 0. 






( III-21 ) 



-vr 



Utilizing the above formula and the convolution 
theorem, recognizing that 



v Q (X 1 ) =0, X 1 < 0 , 



(III-22 ) 



and reintroducing X after inversion, yields 



Q(2np+p-Y) 



Kl(X-S) 

^(X,Y) = -^ / ^ E {v (S)e ia [e~* + 

-1 -VX^ n=0 U 

2 



Q( 2np+ p+Y ) 
1 (X-S) 



Q( 2np+2 p+Y ) ~ Q(2np-Y) 

] - v Q (S)[e * U " S) + e 15 

(HI-23) 



Then, for the potential on the reference blade 
upper surface. 



ifj(X,0 + ) = 1 



X ~(X-S) » 



/ e 



Q 



Q(2np+p)‘ 



- v, 



VttQ' -1 ^ j x _ s ' n=0 

Q(2np+2p)^ Q(2np)^ 

/Site'’ (X - S) + e ' 1 nx^ 7 ])dS . 



Z (v.isie^tse ftsr 3 



(111-24) 



This equation was programmed for digital calculation (see 
Appendix A and Section IV-E). 

3. Modification for Sonic Free Jet Problem 

The sonic free jet is of interest as an extreme case 
of wind tunnel wall porosity. That is, the porosity is so 
great that the wall has vanished and the airfoil is in a 
sonic jet constrained at its boundaries not by a zero 
velocity condition but zero pressure differential. 

The boundary condition on the reference blade remains 
the same as before, equation (111-4). The adjacent blade is 
replaced by the edge of the sonic jet upon which the following 
pressure perturbation condition is imposed 



Transforming and substituting as before, the second 
boundary condition becomes 



Then, using both conditions, the solution to the 
transformed equation is 



iKi/i(X,Y=p) + ^ x (X,Y=p) = 0. 



(III-25) 



(iK + p)iMp,p) = 0. 



(111-26) 




v 0 (p) 



,- sinh[y( Y-p) ] -. 
1 cosh(yp) J 



(III-27) 



y 
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With 



1 = 2e~ yp 

cosh(yp) ^ i+e - 2 p P ) 



2e 



-yp 



" (.l) n e -2MP , 
n= 0 



(111-28) 



*(P,Y) 



V 0 ( P) 

y 



E (_l) n {e" y ^" Y+2np+p ^ 

n=0 



e -y[Y+2np] 



}. (III-29) 



The subsequent inversion is similar to that in the 
previous section. 

The result was programmed for digital computation and 
used as a check case against which the Collocation technique 
porous wall wind tunnel solution was compared. The Fourier 
Transformation method (see next section) was also applied to 
the free jet problem and, by following a derivation similar 
to that above, it was found that an interblade phase angle 
of zero in the cascade problem correctly modeled the sonic 
free jet problem. 

B. SOLUTION THROUGH THE FOURIER TRANSFORMATION 
1 . Cascade Geometry and Boundary Conditions 
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b c n 




Fig. 2. Unstaggered Cascade. 



The arrangement of the cascade is the same as the 
previous section, however, the coordinate system is changed 
to simplify the derivation. The origin of the coordinate 
system is now fixed to the reference blade leading edge. 
Blade chord is c. 

The following nondimensional variables are defined 
for this phase of the study: 



X 





h 




y 



(III-30) 



A 



1 




3 ^ 




All variables will use c (vice b) for nondimension- 
alization as appropriate, i.e., 



K-l = gr . (in-31) 

The boundary conditions are the same as those expressed 
in the previous section, though their representation is 
somewhat different. 



^ y (X,Y=0) = w(X) , 

(X, Y=h) = e ia w(X) , (III-32) 



and 



w(X) = [-l-iK i (X-A 1 )]a 0 + iK ^ . (HI-33) 

2 . Transformation of the Sonic Flow Equation 

As before, begin with the equation for oscillatory 

f 

sonic flow from Section II and transform it with the 
following Fourier Transform, 

-f 00 

f(v,Y) = / f(X,Y)e lVX dX , (111-34) 
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Obtain 



and 



— 2 — 

4>yy + y 'J' = 0 , 



( HI-35 ) 



y 2 = - iXv + 2K x v + K ± 2 



(111-36) 



This equation is satisfied by 



\p(\>,Y) = E,e lyY + F e“ iyY 



(IH-37) 



Using equation (111-32) 



^y(v,h) = e% y (v,0) , 



E 1 [e ilJ ^ - e i0 ] = P 1 [e- il,h - e l0 ] , 



or* 



, -ia iyh 

E i ■ 

1 - e e 



( 111 - 38 ) 



Then, with 



if> v (v,0) = w ( v ) , 
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